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Abstract 

Many classical encoding algorithms of Vector Quantization (VQ) of image compression that 
can obtain global optimal solution have computational complexity O(N). A pure quantum VQ 
encoding algorithm with probability of success near 100% has been proposed , that performs 
operations 45\/iV times approximately. In this paper, a hybrid quantum VQ encoding algorithm 
between classical method and quantum algorithm is presented. The number of its operations is 
less than y/~N for most images, and it is more efficient than the pure quantum algorithm. 
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I. INTRODUCTION 



The two open research problem of vector quantization (VQ) are to design real-time en- 
coding algorithm and to design the codebook generation algorithm that can generate ap- 
proximately local optimal solution (finding global optimal solution is NPC problem) [1]. The 
property of local clustering is revealed by Pang [1] and it is used to design fast codebook 
generation algorithm, that is faster than the famous LBG algorithm by a factor of 4~13 [2]. 
But the other open problem, i.e., fast encoding of VQ with global optimal solution has not 
been solved. The encoding is described as [1]: 

Let codebook 

C = {c[i] | < i < N, c[i] e R\ k > 1} 

For an arbitrary input vector x G R k , we have to find the index i such that distance 
d(x, c[io\) is minimum, where R k is k— dimensional Euclidean Space. Vector c[i] is called 
code vector. 

Grover's algorithm [3] is famous quantum search algorithm, and its accuracy of rotation 
is improved recently [1]. G.-L. Long proposes a modified Grover's algorithm [5]. Long's 
algorithm has the probability of success 100% even for the case that the size of database N 
is very small. This property of Long's algorithm is suitable for quantum image compression 
because N is not a giant number in image compression in general. Some hybrid search 
algorithms have also been proposed to improve efficiency of algorithms [6, 7], in which more 
physical resources are used to commute the efficiency. 

Boyer, Brassard, H0yer and Tap propose an iteration algorithm named BBHT algorithm 
in this paper, in which the Grover iteration is applied repeatedly [8]. BBHT algorithm 
(Boyer, Brassard, H^yer, Tap) [8] is designed specially to solve the searching case that the 
number of solutions is unknown. 

Quantum image compression is possible [9, 10, 11]. Pang et.al., present a Quantum 
Discrete Cosine Transformation (QDCT) for image compression [9]. The method of QDCT 
can be applied to other transformation such as Fourier transformation and classical output 
will be brought. Latorre presents a nice and maybe hopeful quantum image compression 
method. Latorre's algorithm requires that quantization method (scalar quantization or VQ) 
is incorporated with it. Otherwise, it is not competitive with classical methods [10]. Pang, 
et.al., present a pure quantum VQ encoding algorithm with probability of success nearl00%, 
that performs operations 45\/N times approximately [11]. 

This paper presents a hybrid quantum VQ algorithm between classical method and quantum 
algorithm to solve the open problem of fast encoding of VQ. And the number of its operations 
is less than y/N and its probability of success is 100% approximately. 

Notation 1 x: input vector 

Notation 2 S = min{d(c[i], c[j}) \ i ^ j, < i, j < N} 
Notation 3 /: whole encoding space 

Notation 4 S = {x | d(x,c[i ]) < 4f} andT = {x | d(x,c[i )) < 5}, where 5 > 4f. 
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Notation 5 Q(c[z']) = [c\j] \ d(c[i\,c\j\) < 2S,c[j] E C} 
Notation 6 Info = min{| Q[c[z']] || < i < N} 

All of the notations are diagrammatized in figure 2. 

Parameter 5q can be calculated before encoding. The classical data structure of fi(c[i]) 
can be constructed using classical method before encoding (see figure 3). 

II. QUANTUM VQ INEQUALITY ITERATION AND QUANTUM VQ ENCOD- 
ING 

Pang et.al., present a quantum representation of image to which database technique is 
applied [9] , and this quantum method is also used to represent data of image in this paper. 

Definition 7 Oracle Od : 

\5) \x) \i) \c\i\) |0) ^ \6) \x) \i) \c\i])\d(x,c[i\)) 
Od is used to compute the distance d(x, c[i]) and (Od) -1 is the inverse transformation of 

it. 

Definition 8 Oracle Of : 

\S) \x) \i) \c\i}) \d(x,c\i})) (-l)'W \S) \x) \i) \c[i\) \d(x,c\i})) 
, where function f{i) is defined as 

fU\ = { 1 if d{x,c\i]) < 5 
1 ^ ' y otherwise 

Definition 9 quantum VQ inequality iteration Gi nequa ut y : 

G mequality = (2 \o (el - i){u L y\o d y l o f o d u L 

N-l 

, where |£) = \i) and the unitary operation Ul and (t^)" 1 is LOAD operation that 

i=0 

loads data into registers from memory [9]. 
Definition 10 operation Ul: 

\5) \x) |<> |0) |0) U -i \5) \x) \i) \c\i\) |0) 
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A. Sub-procedure 1: the Quantum VQ for x G S 

Algorithm 11 Sub-procedure 1 
Step 1 Let 5 = f . 

N-l 

Step 2 Generate the initial state] ^ >= ~m (1^) \ x ) K) 1°) I ))- 



/N . 
i=0 



Step 3 Perform the iteration G ineqU aUty t 



IV N 



times: 



| \&i >— (G inequality) 1 | IpO > 

Step 4 Observe the third register: let i be the outcome. Compute d = d(x, c[i }) classically. 

If d < 5 , the result is io , otherwise go to Sub-procedure 2. 

If there exists a vector c[i ] e C such that <i(x,c[zo]) < the vector c[i ] is the closest 
codevector of x G 5 and it is the unique closest codevector [1]. And the iteration G inequa ii ty 
acts only on N— dimensional subspace [9, 11]. Therefore, by Grover's algorithm, sub- 



procedure 1 performs the operation Gi nequa ut y 



times approximately. That is, sub- 



procedure 1 has time complexity 



B. Sub-procedure 2: the Quantum VQ for x G (T - S) or x G (I - T) 

Algorithm 12 Sub-Procedure 2 
Step 1 Initialize m = 1 and A — |. 

Step 2 We choose an experiential value 5 such that the set T can cover the whole encoding 
space / approximately. 

Step 3 Choose j uniformly at random among the nonnegative integers not bigger than m. 

N-l 

Step 4 Generate the initial state | $ >= 4= (\$) \ x ) K) 1°) I )) 

i=0 

Step 5 Apply j iterations of the iteration Gi nequa ut y starting from the state | > : 



| ^ j > — (G inequality)' 1 | ^ 

Step 6 Observe the third register: let h be the outcome. 
Step 7 Perform classical computing yo = d(x,c[h]). 

If yo < 5, classically full search the corresponding neighborhood f2(c[/i]) to obtain the 
closest codevector c[io] ; The result is io and exit. 
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Step 8 Otherwise, set m to min{Am, \^N} and go to step 3. 



Step 9 Encoding for I — T : If y > 5 or i does not exist, apply classical method to encode. 



We have conclusion that the closest codevector c[io] of vector x G T — S is included in the 
neighborhood f2(c[/i]) of step 7 (see figure 1). The reason is that d(x, c[j}) > d(c[j],c[h]) — 
d(c[h],x) > 5 for an arbitrary c[j] Q(c[h}). 

The hybrid idea of sub-procedure 2 is that applying the iteration G inequa uty repeatedly 
to find a neighborhood fi(c[/i]) which comprises very few elements, then full search the 
neighborhood classically to find the solution. 

In addition, the set / — T comprises very few vectors statistically and classical method is 
applied to encode vector x £ I — T. 

Figure 1 illustrates the hybrid idea of sub-procedure 2. 



Insert Figure 1 Here 

The iteration G inequa nty is embedded in the BBHT algorithm in sub-procedure 2. The 
case of no solution is handled by BBHT algorithm [S] . The iteration Gi nequa ut y acts only on 
N— dimensional subspace [9, 11]. Thus, the number of performing the iteration Gi nequa uty 

in sub-procedure 2 is not bigger than | \Jj^j^ when Infn <C N. 

The four statistical properties of image for VQ are listed below (see figure 2): 



Property 1 Almost feature vectors x concentrate on themselves centroids generally. That is, almost 
vectors x£l are included in set S. The solution is unique for x G S. 

Property 2 The set T — S comprises the points of near edges of image in practice. The solution 
is included in a small neighborhood f2(c[/i]) for input vector x G T — S . 

Property 3 The set / — T comprises special points such as very detail points or points at edges 
maybe. 

Property 4 Let \ S \: \ T — S \: \ I — T \ = a : b : c. Statistically, we have (see figure 2) 

80% < a < 99%, 1% < b < 19%, and c < 1% 

Therefore, the time complexity of the whole quantum encoding algorithm is less than 
VN statistically. 

The phenomenon should be noticed, that the above properties are not powerful for 
classical methods to solve the open problem of the fast encoding of VQ, by contrast, it is 
powerful for quantum methods. 



Figure 2 shows these statistical properties of image. 



Insert Figure 2 Here 

The hybrid algorithm is not only an image compression method, but also it can be applied 
to image recognition [12]. 



5 



III. THE COST PHYSICAL RESOURCE FOR ACCESSORIAL DATA STRUC- 
TURE 

The main accessorial resource is cost to save the classical data structure of all neighbor- 
hoods f2(c[z]). And the data structure can be constructed using adjacency list technique 
[13] before encoding. Adjacency list technique is a basic technique, and it is applied on 
image compression usually [1] . The adjacency list to save all f2(c[z]) can be constructed by 
the following method [13]: 

First, Allocate a single list to save all neighbors of each codevector c[i], where c[i] € C(i.e., 
centroid in figure 2) and % = 0, 1, (N — 1). 

Second, All of the addresses of all single lists are saved in an array. 
The structure of the adjacency list is diagrammatized in figure 3. 

It is easy to estimate the space complexity of the adjacency list' 13 '. The total number 

N-l 

of the classical bits required to save the adjacency list is approximately (I ^( C KI) I 

+ l)(log 2 iV + 4 x 8) (bits). In general, | fi(c[z]) | is very small (see figure 2). 

Therefore, the space complexity is 0(Nlog 2 N) bits approximately and it is very good. 

It's very easy to realize the adjacency list for modern electronic computer [13]. Therefore, 
it's a good tradeoff that costing 0(N \og 2 N) classical bits to permute the running time by 
decreasing 45 factor than the pure quantum algorithm presented in Ref . [1 1] . 

Insert Figure 3 Here 

IV. CONCLUSION 

VQ has high compression ratio and simple structure. However, the performance comes 
at the cost of increased computational complexity. Fast encoding of VQ is an open problem 
and is one of two key techniques of VQ. Many classical encoding algorithms of VQ that 
can obtain global optimal solution have computational complexity O(N). In this paper, a 
hybrid quantum VQ encoding algorithm between classical method and quantum algorithm 
is presented. The number of its operations is less than y/N for most images. It's a good 
tradeoff that costing 0(N log 2 N) classical bits to permute the running time by decreasing 
45 factor than the pure quantum algorithm in Ref. [11]. 
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V. FIGURES 
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FIG. 1: The illustration of the hybrid idea of sub-procedure 2: 

1. Embed the iteration Gi nequa m y in BBHT algorithm and apply the iteration repeatedly. This 
quantum method will find a neighborhood £l(c[h]). 

2. The global optimal solution c[io] £ 0(c[/i]), and there are very few elements in f2(c[/i]). Therefore, 
we can classically search J)(c[/i]) to obtain solution. 
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FIG. 2: The illustration of the properties of image that is applied to design the hybrid algorithm: 
(Image Lena with size256 x 256 is divided into 2x1 image blocks to form vectors [1].) 

1. Almost feature vectors x concentrate on themselves centroids generally. And the solution is 
unique for x £ S 

2. The set T — S comprises the points of near edges of image and it includes few points in practice. 
And the solution is included in a small set f2(c[7i])for x 6 T — S. 

3. The set I — T comprises special points such as very detail points or points at edges maybe. 

4. Statistically, | S |»| T — S \»>\ I-T\ 

The above four statistical properties of image for VQ are applied to accelerate the quantum algo- 
rithm: 

1. Sub-procedure 1 acts on the set S with time complexity f v~N approximately. 

2. Sub-procedure 2 acts on the set T — S with time complexity less than approximately. 

3. Classical full search algorithm acts on I — T. 

The phenomenon should be noticed, that these properties are not powerful for classical methods 
to solve the open problem of the fast encoding of VQ, by contrast, it is powerful for quantum 
methods. 
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FIG. 3: Fig. 3 The illustration of the adjacency list to save all neighborhoods fi(c[i]): 
Each single list saves all neighbors of c[i}. The space complexity of the adjacency list is 0(N log 2 N) 
and it is very good. It's very easy to realize the adjacency list for modern electronic computer. 
Therefore, it's a good tradeoff that costing 0(A r log 2 N) classical bits to permute the running time 
by decreasing 45 factor than the pure quantum algorithm in Ref.[ll]. 
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This figure "Fig 1. PNG" is available in "PNG" format from: 



http://arXiv.org/ps/cs/0605002v3 



This figure "Fig2.PNG" is available in "PNG" format from: 



http://arXiv.org/ps/cs/0605002v3 



This figure "Fig3.PNG" is available in "PNG" format from: 



http://arXiv.org/ps/cs/0605002v3 



